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Abstract 

We study the relation between the barrier conductance and the Coulomb 
blockade peak splitting for two electrostatically equivalent dots connected by 
tunneling channels with bandwidths much larger than the dot charging ener- 
gies. We note that this problem is equivalent to a well-known single-dot prob- 
lem and present solutions for the relation between peak splitting and barrier 
conductance in both the weak and strong coupling limits. Results are in good 

qualitative agreement with the experimental findings of F. R. Waugh et al. 
PACS: 73.20.Dx,71.45.-d,73.40.Gk 
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I. INTRODUCTION 



Turning on a tunnel junction to a quantum dot leads to progressive destruction of the 
single-dot Coulomb blockade. Experiments by Waugh et al. 0] and Molenkamp et al. || 
chronicle this eradication for two tunnel-coupled dots of equal and widely disparate charg- 
ing energies, respectively Inspired by the experimental results of Ref. the present paper 
seeks to develop a simple model for the coherent tunneling of electrons between a pair of elec- 
trostatically identical quantum dots. (See Figure 1(a) for a schematic view of Waugh et al.'s 
double-dot structure.) The goal is to describe the evolution of the Coulomb blockade from 
that of two isolated dots to that of one composite dot in terms of parameters that deter- 
mine the states of the isolated dots and the nature of the connection between them. In the 
limits relevant to the experimental situation in Ref. 0, we find that the most important 
dimensionless parameters are the number N c h of conducting channels between the two dots 
and the dimensionless interdot barrier conductance g of each channel, which is measured 
when the Coulomb blockade has been removed. (The interdot barrier conductance was mea- 
sured in Ref. by de-energizing the external barrier potentials V x % that separate the dots 
from the leads. This conductance is to be distinguished from the conductance measured in 
the double-dot Coulomb blockade measurements, which will be referred to as the Coulomb 
blockade conductance or double-dot conductance.) 

The problem of coupled quantum dots and more generally, of the effect of tunnel- 
couplings upon the Coulomb blockade has received much attention. I. M. Ruzin et al. Q 
examined the Coulomb blockade structure of two non-identical dots in series via a stan- 



dard activation-energy approach. C. A. Stafford and S. Das Sarma [tl3| , [T5f| as well as 
G. Klimeck et al. JE4J have applied Hubbard-like models with and without interdot ca- 
pacitances to determine the many-body wavefunctions for tunnel-coupling between a small 
array of single-dot eigenstates. Many investigators have studied the effect of tunneling upon 
the Coulomb blockade for metallic junctions, in which there are a large number of con- 
ducting channels [5-9,17]. Relatively few have considered junctions with only one or two 
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channels [16,18-20]. 

In Section II of this paper, we present a brief review of the experimental results that have 
motivated our investigation. In Section III, we define a tunneling model which is useful for 
calculations in the limit of weak coupling between the two dots. The strong-coupling limit 
is analyzed in Section IV, and a summary of our conclusions is presented in Section V. 

II. MOTIVATION 

The experiment of F. R. Waugh et al. [[| provides the primary motivation for this paper. 
These authors study the effect that varying the interdot potential barriers has upon the 
Coulomb blockade conductance peak structure for arrays of n dots, where n equals 2 or 
3. For their Coulomb blockade measurements, they energize the confining gates (V X i in 
Figure 1(a)) so that the conductance between the dots and the external leads is much less 
than 2e 2 /h. Having tuned the dots to be electrostatically identical — i.e., to have common 
gate and total dot capacitances C g and Cs — they find that lowering the interdot barriers 
results in interpolation between the peak structure characteristic of the isolated individual 
dots and that characteristic of a single composite dot having capacitance nC^\ the initial 
isolated-dot peaks split into bunches of n sub-peaks, and the splitting within the bunched 
sub-peaks increases until they are essentially equally distributed with n-times the periodicity 
of the original peaks. (See Figures 2(b) and 2(c).) For the double dot {n = 2), Waugh et al. 
also measure the conductance Gb of the barrier between the two dots after the exterior 
walls of the double dot have been removed. In plots of the sub-peak splitting and barrier 
conductance as functions of the barrier gate voltage, they remark that the latter appears to 
be a displaced repetition of the former, the displacement being explicable as a result of the 
need to correct for the exterior walls' influence upon the barrier. 

Waugh et al. use a T = "capacitive charging model" to interpret their data. In this 
model, electrons on the dots are treated as charged particles with no kinetic energy that 
occupy each dot in integer amounts. In the absence of coupling, the energy is given by the 
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sum of the potential energies of the individual dots. For two dots with common capacitances 
Cs and C g , the expression for the energy has a particularly simple form: 

i=l 

where U is the charging energy for each individual dot, U = e 2 /C*s; rti is the number of 
electrons on the ith dot; and fa is the gate voltage parameter that determines the energy- 
minimizing value of rii. For common gate voltages and gate-to-dot capacitances, we have 
the relations fa = C g iV gi /e = C g V g /e = (f>. Figure 1(a) should help put these parameters in 
context. 

For each set of integer occupation numbers (711,712), the capacitive charging model with 
fa = <p gives an energy Er ni}Jl2 \ that is a parabolic function of the common gate voltage 
parameter 0. (See Figure 2.) All the parabolas are identical in shape, their only distinguish- 
ing features being the locations of their minima. The lowest-energy parabola E^ tot {fa) for a 
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given value of N tot = n, has ri\ = n 2 = N tot /2 for N tot even and ri\ = n 2 ± 1 = Nto ^ ±l 

i=l 

for Ntot odd. In the former even case, the minima all lie on the line E = 0. In the latter 
odd case, the minima are displaced upward, sitting along E = U/4. For all parabolas, the 
0-coordinate of the minimum is N tot /2. 

A prominent peak in the double-dot conductance occurs at values of <p such that the 
lowest-energy parabolas corresponding to consecutive values of N tot cross — in other words, 
at values of (p for which E^ tot {fa) = E Ntot+ i(fa) for some integer Ntot. F° r the model of 
Equation ([]]), this occurs whenever <fi = m + |, where m is an integer. (One such crossing 
point is marked by the black dot in Figure 2.) 

In a model in which coupling between the dots is included, the lowest-energy parabolas 
for odd N to t are shifted downward relative to the lowest-energy even-N to t parabolas by an 
"interaction energy" E int . This downward shift splits each of the initial crossing points 
into a pair of crossing points symmetric about the position of the initial degeneracy, from 
which they are separated by a distance proportional to E int . As a result, each of the initial 
conductance peaks is similarly split into two sub-peaks with separation proportional to Ei nt . 
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The sub-peak splitting reaches its saturation value when E{ nt = U/4 — i.e., when the lowest- 
energy even and odd parabolas have the same minimum energy. At this point, the relevant 
crossing points occur for = \(m-\- 1). The corresponding conductance peaks are once again 
equally spaced, but their period is now that characteristic of a single dot with capacitance 
2C E . 

Thus, in the capacitive charging model, the problem of explaining the peak splitting 
reduces to the problem of describing the shift in the ground state energy of a double dot 
containing a fixed total number of particles. Waugh has shown that introduction of a ca- 
pacitive coupling Ci n t between the two dots would allow one to obtain a picture in qualitative 
agreement with the experimental results: as the interdot capacitance goes to infinity, Ei nt 
converges to U/4. However, the magnitude of the interdot coupling necessary to fit the ex- 
perimental data is much larger than what one would expect from an electrostatic interaction 
between two adjacent dots having a narrow tunneling channel between them. Consequently, 
the use of the interdot capacitance Ci nt must be regarded as simply a reparametrization 
of the problem which replaces one unknown, Ei nt , with another unknown, Ci n t. What we 
really want is a theory which produces at least qualitative agreement with experiment and 
expresses E int in terms of simple measurable quantities. Waugh et al. provide one candidate: 
the conductance Gb of the barrier between the two dots. The remainder of this paper is 
devoted to developing a theory of the relation between E int and the dimensionless conduc- 
tance per tunneling channel g = AT b n , where N c h is the number of independent interdot 
tunneling channels (assumed to have identical conductances) and Go is the conductance 
quantum e 2 /h. In the experiment of Ref. 0, there is no applied magnetic field and the dots 
are connected by a narrow constriction allowing only a single transverse orbital mode with 
double spin degeneracy. As a result, in this experimental case, N ch equals 2. 



III. TUNNELING MODEL FOR THE DOUBLE-DOT COUPLING 



5 



A. Definition of the Model 



Our goal can be stated a bit more precisely. For a general tunnel-coupling between two 
dots involving any number N c h of identical, independent channels and dimensionless channel 
conductance g, our aim is to express the fractional energy shift / = jfjf as a function of g 
and N c h plus any other parameters that might be found to be important. In order to derive 
an equation for /, we first choose a double-dot Hamiltonian. Since we hope to explain the 
evolution of the double-dot Coulomb blockade through the barrier conductance alone, we 
will ignore — for the moment at least — electrostatic coupling of the dots. Interaction between 
the dots will occur solely via tunneling through the barrier between them. Such tunnel 
Hamiltonians have been found useful from the beginnings of Coulomb blockade theory ||11|| , 
and the model we will use is a double-dot version of the Hamiltonian used, for example, 
by Averin and Likharev to investigate the conductance oscillations of small metal-to-metal 



tunnel junctions 12 . In particular, we have the Hamiltonian H = H + Hf, where 



H = K + V, 

2 

k = 2J ^2 ^2 e *°-™*°-> 

i=l a k 
i=l 

H T = ^2Y1 (*kik a 4k a «r c lki«T + h - c -)- ( 2 ) 

<t kik2 

In these equations, i is the dot index, a is the channel index (which could signify different 
spin channels), and k is the index for all internal degrees of freedom not included in the 
channel index. In addition, hi = yjrWk<x ls ^ ne number operator for the zth dot, and tkik 2 

ko- 

is the tunneling matrix element between a dot 1 wavefunction indexed by ki and the dot 
2 wavefunction lying in the same channel and indexed by k2. The gate voltage parameter 
(f>i has the same meaning as in Equation ([!]). is the kinetic energy of the single-particle 
eigenstate of the ith dot having the indicated degrees of freedom. For simplicity, we will 
take these energies to be independent of dot and channel: e^ a = e^. 
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The next step in focusing upon a model Hamiltonian is to choose a form for tkik 2 - 
Quite generally, tkik 2 wm De nonzero only when both k\ and ki lie within some wavevector 
shell that maximally spans the space between the theory's low and high momentum cutoffs. 
The size of the wavevector shell depends on details of the barrier: for a channel with an 
abrupt delta-function barrier, the shell spans the largest possible energy range; for a slowly 
developing adiabatic barrier, the shell width will be small on the scale of a Fermi wavevector. 
Important questions are how many states lie within this shell — i.e., how large is the width W 
of the corresponding energy shell compared to the average level spacing 5 between different 
states in the same channel (hereafter referred to as "the average level spacing" or just "the 
level spacing") — and for a given ki, for how many k 2 is tkik 2 nonzero. Thin-shell models 
with "one-to-one" hopping elements (i.e., for which tkik 2 = unless k x = k 2 ) have been 
applied to the coupled dot problem with some success [|13| HT5|1, particulary for level spacings 



S which are on the order of the charging energy U. For the nearly micron-sized dots used 
by Waugh et al. |2|, however, U is approximately 150 /xeV and 5 is on the order of 10 fjteV, 
so we expect that a tunnel-coupling sufficient to destroy the isolated-dot Coulomb blockade 
will involve a large number of single-dot eigenstates. 

Consequently, we consider here a thick-shell model which is the antithesis of the injective 
thin-shell model. Working in a regime where W ^> U ^> 5, we use a tunneling matrix element 
t that is independent of k\ and k 2 within the shell: 



''ki^ 



t : V ki, k 2 s.t. e < e kl , e ka < e + W. 



As the quantities we calculate are independent of the phase of t, we guiltlessly choose t to 
be real. This model is roughly equivalent to one in which each dot is represented by a tight- 
binding lattice with intersite hopping elements of order W/5 and where interdot tunneling 
occurs via a tunneling Hamiltonian with a single site-to-site connection. Choosing these 
tunneling sites to be at the origins Oi and 2 of the respective lattices, we may write 



h t = ^2(Tcl 02(7 c 10ia + h.c), 



where T = Ny/t and Nw = W/5 is the number of orbital states per channel in each dot 
within the bandwidth W . (The equivalent lattice model should include second and further 
neighbor hopping so that the density of states is approximately constant between eo and 
eo + W. The lattice constant is chosen by requiring that the product of Nw and the area of 
a unit cell equals the area of a single dot.) 

As the Fermi energy ep must be somewhere between eo and eo + W, the meaning of eo 
depends on the width of the band. For a maximally thick shell, eo lies at the bottom of 
the conduction band, and W is an ultraviolet cutoff, chosen to be of order twice the Fermi 
energy. Alternatively, when the barrier between the dots has a broader spatial extent, the 
energy shell sits more narrowly about the Fermi energy, and the width W is on the order 
of the energy difference needed to produce a factor-of-two change in the magnitude of the 
transmission amplitude for an incident particle. We define a dimensionless filling parameter 

_ e F -e 
F = _ ^' 

which gives the position of the Fermi level within the bandwidth W. Provided that (1 — F)W 
and FW are both large compared to U, our final results should be independent of the precise 
values of W or F. 

The model we have constructed is basically the two-dot version of that used by 



L. I. Glazman and K. A. Matveev jn| and H. Grabert |17j to study the charge fluctua- 
tions of a single metal particle connected via point-tunnel junctions to conducting leads. 
(See Figure 1(b).) Indeed, the similarity between the two-dot and one-dot problems is even 
more fundamental than this observation indicates. Consider again the double-dot potential 
energy V. By transforming to the analog of center-of-mass coordinates, one generates the 
following form: 

V = ^(N tot - $ tot ) 2 + U{n - p/2) 2 , (3) 

2 2 

where N to t = yj^i> &tot = n — n2 l 2 ni , p = 02 — 01- The rationale for the nor- 

i=\ i=l 

malizations for n and p will soon be made apparent. In the meantime, note that for our 



Hamiltonian, N tot is a constant of motion. Thus, for given N tot , &tot, and U, we can drop 
the first term and insert in the Hamiltonian a reduced potential energy: 

V red = U{n-p/2f. (4) 

Now restrict N tot to be even. Then, n has integer expectation values in all the un- 
perturbed double-dot eigenstates. The Hamiltonian is exactly that of a single dot tunnel- 
coupled to an ideal lead. The dot has number operator h, charging energy 2U, and gate 
voltage parameter p/2. In the absence of tunneling and with the level spacing in both dots 
much less than the charging energy U, the ground state is an eigenstate of n that minimizes 
the reduced potential energy, which in the future we consider equivalent to "the potential 
energy." For p = 0, the minimum potential energy is zero and is achieved when the eigen- 
value n of n is zero — i.e., when there are an equal number of particles in the two dots. All 
other values of n give higher potential energies. For p = 1, on the other hand, the minimum 
potential energy is U/4, and n = and n = 1 give degenerate minima. 

These no-tunneling distinctions between zero and U/A and between nondegeneracy and 
double degeneracy are quite familiar: they characterized the previously discussed even and 
odd double-dot ground states (p = for both). Indeed, what we called the 11 even double- 
dot ground state" is precisely the u N tot even, p = ground state." The "odd double-dot 
ground state" is not exactly the same as the u N tot even, p = 1 ground state"; there is 
no getting around the fact that one case has one more (or less) particle than the other. 
However, in terms of their ground state energies, the difference between the two will be 
down by a factor of FN W or (1 — F)N W . For a wide shell somewhere in the vicinity of 
half-filling, both FN W and (1 — F)N W are much greater than one, and the above difference 
is negligible. Calculation of E int with 4>i — 02 for the double dot is therefore effectively 
equivalent to calculating the relative shifts of the p = and p — 1 ground states of a single 
dot tunnel-coupled to a bulk two-dimensional electron gas. More generally, the value of p 
can be arbitrary, and the difference in the ground-state energies of the double-dot system 
for even N tot and odd N tot is related to the difference in the ground-state energies of the 

9 



single-dot system for gate voltage parameters p and 1 + p. 

We have mapped our double-dot problem onto a more general single-dot problem for 
which we can calculate the relative downward shift of the p ^ ground state to the p = 
ground state. As the minimum potential energy is periodic in p with period two and is 
also even in p, we need only calculate up to p = 1. Dividing by the zero-tunneling energy 
difference of the two ground states, we find that our emended aim is to calculate 

where < p < 1, u — U/W, Nw = W/5, and E int (p) is the ground-state energy relative to 
the ground-state energy for p = 0. 



B. Barrier Conductance in the Weak-Coupling Limit 

Before we can derive our equation for f p in terms of g, we must find a formula for the bar- 
rier conductance. Measurement of the barrier conductance Gb with the exterior gates turned 
off can be modeled by calculating the tunnel junction conductance for U = 0. As mentioned 
before, we assume the different conducting channels to be identical yet independent — their 
individual conductances are the same and they do not interfere with one another. These 
assumptions are certainly reasonable for the two spin channels in the experiment of Ref. || . 
Using the Lippmann-Schwinger equation with inserted for the scattering potential UlO 



one can solve for the perturbed electron eigenfunctions. The Heisenberg equation of motion 
for ni can then be used to solve for the particle flow from dot 1 to dot 2 for a given voltage 
bias. Solving the resulting expression for the linear conductance gives the following equation 
for the dimensionless conductance per channel: 

9 N ch G ~ |l + x«r 1 ' 

where a = (ttT/W) 2 = ijit/5) 2 and x = (1 + ^og(j^)) 2 . H. O. Frota and K. Flensberg 
have derived this equation for half-filling (F = 0.5, x — 1) v i & a Green's function-Kubo 
formula approach [jlSfl . 
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The calculated conductance Gb exhibits rather curious behavior: it first rises to a maxi- 
mum of N c hGo corresponding to N c h fully open channels (It does this even when im[x] 7^ 0.), 
and then falls asymptotically to zero as (T/W = t/8) — > 00. As Frota and Flensberg 



note |18[], the asymptotic damping of the conductance results from the fact that formation 
of bonding and anti-bonding states at the tunnel junction makes the cost of passing through 
prohibitively high. The limit of (T/W = t/8) — > 00 is in some sense unphysical: we do 
not expect a point-to-point hopping coefficient T to significantly exceed the tunneling shell 
width; nor do we expect the tunneling matrix element t to be much greater than the average 
level spacing. Nevertheless, the apparent absence of any good reason to truncate the theory 
at a particular value of t leaves us with a model that is at best unwieldy in the limit of 
strong coupling. To get the correct limiting behavior for strong coupling, it is more conve- 
nient to use a different approach, suitable for perturbation about the g = 1 limit. This will 
be described in Section IV. 



C. Relative Energy Shift of Even and Odd States in the Weak-Coupling Limit 

In the meantime, the site-to-site tunneling model is still useful in the weak coupling 
regime. So we plod ahead, calculating via standard Rayleigh-Schrodinger perturbation the- 
ory the second order shift in the ground state energy for p 7^ minus that for p = 0. The 
p = 1 shift will be taken to equal the limit of the general < p < 1 shift as p —>■ 1 . It might 
be objected — correctly — that this limit fails properly to account for the degeneracy of the 
ground state at p = 1. Such a failing is pardonable, however, for the contributions that are 
left out are all smaller by a factor of FNy/ or (1 — F)N\y from those which are retained. 
Since we assume that t/8 is finite, F is of order |, and Nw is large, the omitted terms are 
negligible. 

For Nyy 3> 1, the perturbation theory sums can be approximated as integrals. Observing 
that u = <S 1, we divide the difference between the second order shifts by Up 2 / 4 to get 
the leading approximation to f p : 
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f? = 4iV ch ^[(l - p) log(l - p) + (1 + p) log(l + p) + 0(V)]/P 2 - (7) 
The second-order term indicates a significant feature of f p : it is even in p. This property has 



been noted by H. Grabert |L7j and results from the fact that at any order of perturbation 
theory, every tunneling process contributing to the energy shift has a twin with the roles of 
dots 1 and 2 interchanged. In any intermediate virtual state with eigenvalue n for n, the 
potential energy is greater than that for the unperturbed ground state by SV(p) = Un(n—p). 
Therefore, when dots 1 and 2 are interchanged, SV(p) — > 5V(—p) for all the intermediate 
states. If we represent one of the twin terms by A(p), the other is A(— p), and we see that 
f p is constructed of sums that are even in p. 

Using the second-order (in t/S) parts of g and f p , we can now write a first-order equation 
for fp in terms of g: 

, ( i) N ch g [(1 - p) log(l - p) + (1 + p) log(l + p) + Q(up 2 )} 

fp ~—2 -2 > (8) 

a result consistent with the large- N c h calculation of the effective capacitance of a single dot 
at p = 0. Setting p = 1 to calculate the relative shifts of the original even and odd states, 
we find 



f {1) = ^^N ch g + 0(ug,g 2 ), (9) 



2 log 2 



7T 2 



where we have used the fact that / as originally defined without the subscript is equivalent 
in our limits to f p =\. The above equation indicates that the suggestion that the plot of / 
is a displaced version of the plot for g is not precisely correct. In particular, for g <C 1 and 
N c h = 2, Equation @ gives a slope of approximately 0.28 for /(<?), rather than unity. Thus, 
in this regime, the fractional splitting / of the double-dot conductance peaks should lag g, 
the dimensionless barrier conductance per channel. 



IV. CONNECTION TO THE STRONG-COUPLING LIMIT 



If we blithely extended our perturbative equation for / to the limit g — > 1, the large- N c h 
f would greatly overshoot its mark and the one or two-channel / would fall substantially 
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short. The real issue is not, however, how badly such a naive extrapolation fails, but whether 
we can interpolate between these weak-tunneling results and those which can be calculated 
for the strong-tunneling limit. For the large- N C h limit, a reasonable interpolation between 
the solutions for weak and strong coupling has already been found |)]|)||§. The situation 
is less clear when N c h equals one or two. Flensberg and Matveev |B5],^D| have proposed a 



useful Luttinger-liquid approach in which the nearly transparent link between a dot and 
an electrode is modeled as a one- dimensional channel with a slightly reflective potential 
barrier. Convergence to the single composite-dot limit is achieved naturally and neatly, and 
Ei n t is calculated perturbatively in r, where r is the reflection amplitude, and g = 1 — |r| 2 . 
Translating Matveev's calculations of the leading term for (1 — g) <C 1 into our language, 
we find that for N ch = 1 (i.e., assuming spin polarization), 

f=l--^y/T=]j, (10) 

where 7 ~ 0.577 is the Euler-Mascheroni constant. For the case relevant to the experiment 
of Ref. 0, N c h = 2, Matveev's calculation gives 

f = l+™£.(l-g)]og{l-g). (11) 

Except for the logarithmic factor in the second formula, these equations are of the form 



suggested by the scaling analysis of Flensberg [ 19 1 , which predicts effective charging energies 



behaving as (1— g) Nch ^ 2 . Matveev's initial two-channel solution is, in fact, linear in (1 — (?) but 
diverges logarithmically as U/8 — ■> 00. A higher-order analysis to eliminate the divergence 
replaces the logarithm having argument U/8 with one having argument (1 — g)^ 1 [f20 . 



In Figure 3, we show the /-versus-g plots given by the weak and strong coupling formulas 
(|), (0), and fln|) for N ch = 1 and N ch = 2. In each plausible interpolation between 

the weak and strong coupling limits is given by a dashed curve. The experimental data of 
Ref. are in reasonable agreement with the dashed curve for N c h = 2. Nevertheless, from 
a theoretical standpoint, it is clear that, unlike the calculations for N c h ^> 1, for N C h = 2 the 
order of calculation completed so far does not really allow confident interpolation between 
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the weak and strong coupling limits. Calculation of higher orders in perturbation theory 
should improve the matching, but such computations are made difficult by the fact that 
the correlations induced by the strong Coulomb interaction make normal Green's functions 
methods inapplicable ||. Different time orders must be treated separately, and as appears 
to occur quite generally in Coulomb-blockade problems [21], the number of diagrams grows 
pathologically with the order in perturbation theory. Nevertheless, calculation of the g 2 -term 
in the weak-tunneling limit is within reach, and this term may suffice to give a more reliable 
interpolation between the weak and strong coupling regimes. 



V. CONCLUSION 

Following the work of Waugh et al. ||, we have investigated the relation between the 
barrier conductance and the Coulomb blockade for two electrostatically equivalent dots 
connected by one or more identical tunneling channels. We propose to write the fractional 
peak splitting / of the Coulomb blockade conductance peaks as a function of the number of 
channels N c h and the dimensionless barrier conductance per channel g, assuming that the 
energy level spacing 5 is small compared to the Coulomb blockade energy U and that U is 
small compared to the bandwidth W of states over which the amplitudes for transmission 
through the barrier are roughly constant. Using a "uniform thick-shell model" for the 
tunneling term in the Hamiltonian, we solve for this function to leading order in the limit 
of weak interdot coupling. We find that the peak splitting should evolve substantially more 
slowly than the barrier conductance in this limit. Having noted that the two-dot problem 
can be mapped onto a better known one-dot problem, we adapt previous strong-coupling 
results to obtain the asymptotic form of the double-dot peak-splitting in the limit g — > 1. 
For the case of N c h = 2, which is pertinent to the experimental results of Ref. ||, the 
limiting forms for the strong and weak coupling do not match up well enough to allow a 
reliable quantitative interpolation between the two limits. Nevertheless, a visually plausible 
interpolating curve is in good qualitative agreement with existing experimental data. 
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FIGURES 

FIG. 1. (a) Schematic diagram for the double dot. Negative potentials are applied to each of 
the gates to form the double-dot structure. The gate potentials V g \ and V g 2 control the average 
numbers of electrons on the dots. These are the potentials that are varied to see the Coulomb 
blockade. V& controls the rate of tunneling between the dots. V x \ and 142 control the rate of 
tunneling to the adjacent bulk 2D electron gas (2DEG) leads. For calculations of the double-dot 
energy shifts, tunneling to the leads is assumed negligible compared to tunneling between the two 
dots. In measuring the barrier conductance Gb, however, the potentials V x i are turned off so that 
each dot is strongly connected to its lead. The side-wall potentials V^i and V S 2 are fixed, (b) 
Schematic diagram for the single dot. V& now controls tunneling between the dot and the bulk 
2DEG. V g determines the average number of electrons on the dot. For our purposes, V s and V x are 
constant, and tunneling to the bulk 2DEG through the barrier defined by V x is negligible compared 
to tunneling through the barrier defined by V&. 
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FIG. 2. (a) Energy curves in the capacitive charging model for electrostatically identical dots 
with V g i = V g 2 ■ Energies are given in units of the charging energy U ; the gate voltage is given in 
units of £T- . Each zero-coupling eigenstate with definite particle number n{ on the ith dot gives rise 
to a parabola, labeled (n\,n2), which shows the state's energy as a function of the gate voltage. 
The zero of energy is chosen to coincide with the lowest energy possible for states with an even value 
for the total number of particles N tot . The solid odd-iV toi parabola gives the lowest-energy curve 
only when there is no interdot coupling. The dotted parabola is the shifted-down energy curve for 
odd N to t that results from finite coupling between the dots. The relevant degeneracy points are 
indicated by a black dot for zero coupling and white dots for finite coupling, (b) "Zero-coupling" 
conductance through the double dot as a function of the gate voltage. (For ease of viewing, peaks 
are depicted as symmetric with uniform finite widths and heights.) Conductance peaks are aligned 
with the zero-coupling degeneracy points such as the one shown in (a) and occur regularly with 
unit period, (c) Conductance through the dot for finite interdot coupling. Conductance peaks are 
aligned with the perturbed degeneracy points. Each zero-coupling peak has split into two separate 
peaks, equally distant from the zero-coupling peak position. Increasing the interdot coupling 
increases the separation between the paired peaks until the full set of peaks is again regularly 
distributed, with half the original period. (This figure for the capacitive charging model follows 
that of Ref. [2].) 

FIG. 3. Graphs of the fractional Coulomb blockade conductance peak splitting / as a function of 
the dimensionless conductance per channel g in the weak and strong tunneling limits for (a) N c h = 1 
and (b) N c h = 2. Possible interpolating functions are shown by dashed curves. 
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